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TECHNICAL MEMORANDUM X- <;474‘) 


A SIMPLE PERFORMANCE CALCULATION METHOD FOR LH 2 /LOX 
ENGINES WITH DIFFERENT POWER CYCLES 

SUMMARY 


The selection of a feed system for an engine requires the estimation of 
its effects upon performance. For an engine with a gas generator cycle, the 
performance depends upon the amount of gas generator flow rate for the tur- 
bine drive. Based upon a power balance, the necessary mass flow and pres- 
sure formulations are derived. This leads to a simple equation for the per- 
formance, gas generator mass flow rate, and pump discharge pressures. The 
combination liquid hydrogen/liquid oxygen (LH 2 /LOX) is selected as a typical 
propellant combination for vacuum space flight applications. Approximate 
equations for the performance of the combustion products are derived. The 
solution method has been applied to five different engine types leading to an 
optimization of the chamber pressure and turbine pressure ratio associated 
with maximum performance. 


INTRODUCTION 


The propulsion system has a dominant place in astronautics, since its 
performance essentially decides the feasibility and success of a mission. To 
increase the performance of an engine at sea level, the chamber pressure must 
be increased. This reduces the effects of dissociation, but much more impor- 
tant is the use of a larger area ratio. 

For pure- vacuum space flight the situation is different, since the expan- 
sion ratio is limited only by the available geometric envelope of the engine. An 
engine using a gas generator cycle needs a certain amount of the whole mass 
flow to drive the turbopumps. Although the specific impulse of this type is 
lower than that of an engine with a staged combustion cycle, which uses the 
whole mass flow through the turbine in the main chamber, its use may lead to 
some advantages. The pump pressure for the gas generator system is lower 
and the development costs are smaller. The decision, which type of engine 
should be used, depends greatly on the amount of propellant flow necessary for 



the turbine drive. The performance of the engine with the gas generator cycle 
is very much influenced by the turbine and pump technology. A simple method 
for the calculation of the specific impulse, the necessary gas generator flow, 
and the pump, gas generator and chamber pressure permits an optimization of 
the principal design and a quick estimation of the effects of the various param- 
eters on engine performance. 


BASIC EQUATIONS FOR PERFORMANCE CALCULATION 
Gas Generator Engine 

A schematic representation of a rocket engine with a gas generator 

cycle and the assigned nomenclature is shown in Figure 1. The mass flow of 

oxidizer and fuel at pump pressure (p ) into the gas generator is denoted by 

P 

m and m. , respectively. After the combustion at pressure p the 
ox fu g 

g g 

gases expand in the turbine to a pressure p and are injected into the main 

te 

chamber nozzle. The injection occurs at a position where the main chamber 
nozzle's static wall pressure and the turbine exit pressure are almost the 
same. To simplify the calculations, no mixing and secondary chemical reac- 
tion of the turbine exhaust gases with the combustion chamber gases are 
assumed. The expansion of both mass flows is calculated independently. 

Mass Flow Through the Gas Generator. For the calculation of the gas 
generator mass flow, an energy balance is performed. The required energy 
for the pressure increase of the propellants is 


-P: 


E = rh 
p ox 


ox 


ox 


m 


ox 


fu 


ox 


-P. 

^fu Jfu 

p fu 


( 1 ) 


fu 


In this relation m and m„ denote the whole oxidizer and fuel mass flow, 
ox fu 

p and p. the pump outlet and inlet pressures, p the propellant densities, 

p i 

and 7}^ the pump efficiencies. Assuming for convenience that both pump 
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Figure 1. Gas generator cycle engine. 


efficiencies are equal and that the pump inlet pressures are small compared 
with the discharge pressure, equation (l) can be rewritten as 


E 

P 
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ox 


ox p 


+ m 


ox 


fu 


fu p 


fu 


( 2 ) 
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Introduction of the overall mixture ratio r , 


m 


r = 


ox 


m 


( 3 ) 


fu 


results in 


m 


ox 


p(l + l/r) \ P ox 


fu 1 
r 


fu 


( 4 ) 


The turbine shaft energy (E^) results from the enthalpy drop (Ah^) 
between the gas generator and the turbine exit, 


E t ■= A g Ah t , t 


( 5 ) 


where m represents the gas generator mass flow and r? the turbine effici- 
g t 

ency. Instead of the enthalpy drop, the de St. Venant-Wantzei equation can be 

used* and one obtains 



m tj (c T) 
g t p g 



( 6 ) 


Here (c T) describes the total enthalpy in the gas generator, p the gas 
P g g 

generator pressure, which is equal to the turbine inlet pressure, and p the 


* Since the gas generator mixture ratio is normally very low, dissociation 

effects can be neglected and therefore equation (6) is a reasonable equation 

with constant y . 
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turbine exhaust pressure. The 

the turbine is denoted by y . 

g 


isentropic exponent oi' the combustion gases 
Under steady-state conditions both energies 


in 


must be equal. 


E 

P 


E 


t 


( 7 ) 


This results in a gas generator flow 


m 

S 



( 8 ) 


Replacing the absolute value (m ) by the ratio of the gas generator flow to the 

g 

total flow, 


k 

m 


m 

(f 

& 

m 


0) 


the relative gas generator mass flow in terms of the overall mixture ratio is 


k 

m 




(c T) 

v p 'g 




y 

cr 
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( 10 ) 
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Normally the gas generator mixture ratio (r ) and the combustion 

S 

chamber mixture ratio (r ) are given and r can be expressed by these two. 
With c 
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ox 
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fu 


r = r 
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1 + 
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ox 
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fu 
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e transformations and the 
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1 + l/r 

m 

g 
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1 + r 

m 

c 

1 - k 

1 + r 

m 
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( 11 ) 


( 12 ) 


Due to pressure drop in the lines, injector, etc, , the pump outlet pres- 
sure must be higher than the chamber pressure (p ). A simple relation is 
formulated by 



p c (l 


k ) 
ox 


( 13 ) 



p c (i + 



( 14 ) 


where the k values are assumed to be constant within certain pressure limits. 
Since the results are not very sensitive to the pressure factors, this statement 
contains little error. 
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As gas generator pressure, the oxidizer pump discharge pressure, 
which is commonly lower than the fuel pressure, can be used following the 
previous formulation, 



( 1 - k J 


ox 


- ( 15 ) 


= p(l+k )(1 - k ) 
c ox g 


(16) 


With the knowledge of the different constants for a given chamber pressure, 
the gas generator mass flow with turbine exit pressure as a parameter can be 
calculated. 


. From equation (10) a minimum gas generator flow results. The lowest turbine 
flow exit pressure for injection into the main nozzle is limited by the nozzle 
exit pressure (p^). The minimum gas generator mass flow^k^ ^is deter- 
mined by 


i m 

\ min 


! 


k 

m 

nun 




(c T) 

P g 



y - 

o* 

-J71 




(17) 


The numerical value results from a solution of equations (12) and (17). An 
approximate value can be obtained by using r^ instead of r in (17), since 

the gas generator mass flow is normally small compared with the total flow. 

An increase of the gas generator mass flow permits for the same turbine power 
generation the pressure ratio in the turbine to decrease. 


Engine Specific Impulse. The engine specific impulse (1 \ 

V SP ggc/ 

ing from the main chamber thrust and the turbine exhaust thrust can be 
expressed by 


result- 


7 



( 18 ) 


sp. 


ggc 


(1 - k )I 
m sp 


k I 
m sp 


g 


sp 


where I and I 

sp sp 

c g 

impulse, respectively. 


describe the main chamber and gas generator specific 
r]^ considers the various losses such as divergence, 


sp 

friction, etc. For a given gas generator flow, the turbine pressure ratio can 
be obtained by solving equation ( 10) . This results in 



1 



1 1 + k fu\ 
r ' p fu / 


y -i 

g 

' y„ 

g 


( 19 ) 


The pressure ratio of the turbine gases for the expansion in the main 
nozzle is 






k )(1- 
ox 


k ) 
g 


( 21 ) 


where the main nozzle pressure ratio (p /p ) is defined by the engine expan- 
sion ratio. 
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Staged Combustion Cycle Engine 

In a staged combustion cycle engine, which is schematically represented 
in Figure 2, the turbine exhaust gases are fed into the main chamber through 
the injector plate. In this type of engine the expression preburner, instead of 
gas generator, is used. The calculation of the required turbine pressure ratio 
can be done with the previous equations. For a low pressure drop in the tur- 
bine, according to equation (19), the maximum possible flow rate in the 
preburner has to be used. This means that the whole fuel mass flow is injected 



Figure 2. Staged combustion cycle engine. 
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into the prebumer. With the subscript pb for the preburner, the mass flow 
is 

m, = m c . (22) 

fu , fu 

pb c 


The preburner mixture ratio is defined by 


m 


ox 


J2k 


pb m 


fu 


pb 


and one obtains 


m 

-J& 

m 


1+ r 


Pb 


1 + r 


= k 


m 


(23) 


Then with equation (19 ) the pressure drop in the turbine is calculated. 
Analogous to the linear relations for the pressure drops, a similar equation 
can be set up for the connection between turbine exit and main chamber 
pressure, 


te 


p o ( 1 + 


k J 

pb 


(24) 


The pressure drop in the turbine of a staged combustion cycle is normally much 
lower than that of a gas generator cycle. Therefore, equation (19) can be sim- 
plified by expanding into a power series and using only the linear term, 
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1 + k • 1 + k 

ox , 1 fu 

+ — I r 

P r p / c 

ox c f 


fu 
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pb 


Y 1 (1 + r )(c T) 
t p pb p pb 


(25) 


In equation (25) the preburner pressure has to be used instead of p^ , since 

the pumps supply propellants to the preburner. Then together with equation 
(24), equation (25) can be rearranged and one obtains 
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Pb 

te 
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. p (l + k ) 
7,-1 c pb 
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1 + k 1 + k 

ox 1 fu 

+ ir 
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P.P 0- + r )(c T) 

t p • pb p pb 


The pump discharge pressures are 



p c u + k 




P 


te 


(1 + k ) 


(27) 


and 



/ 


P (1 + k ) (1 + k ) 

c pb p^ fu 

te 


(28) 


In a staged combustion cycle the enthalp}' drop in the turbine is recov- 
ered by the pump head. Therefore, the specific impulse can be calculated in 
the usual way with the tank enthalpies of the propellants. With the subscript 
see for staged combustion cycle engine, the specific impulse is 
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(29) 


sp 


see 


sp. 


sp 


The various losses are comprised in 



Approximate Equations for the Thermodynamic Properties 
of the Combustion Products of O2/H2 

For a simple treatment, a specific impulse relationship using only area 
ratio, chamber pressure, and mixture ratio is convenient. Since exact values 
for one-dimensional expansion can be obtained by Reference 1, it is practical 
to use "approximations to solutions" instead of "solutions of approximate 
equations" [2j. The de St. Venant-Wantzel equation belongs to the latter case 
and has the disadvantage that mixture ratio effects, etc. , are omitted almost 
completely. This leads, especially in the near of the impulse optimum, to 
great discrepancies between these approximate values and the exact ones. 

Therefore, simple equations should be used, which represent the exact 
one-dimensional values within small errors. 

Combustion Chamber Gases. The vacuum specific impulse I of 

SP o 

an engine supplied with LOX/LH 2 at tank enthalpy level for different area ratios 
and chamber pressures is presented in Figure 3. Instead of a pure I 

SP o 

approximation, a combination of the characteristic velocity (c*) and the thrust 


coefficient (C ) is convenient. 

J. 

With C n 

as the conventional thrust coeffici- 


X 

0 


ent the usual relation yields 



* / A 



b P - r is + t 

o o Y o t 

k ) 

(30) 


The use of c* , C^, , and p /p approximations has the advantage that 
r c e 

o 

individual effects can be analyzed. 
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490 



Figure 3. Vacuum specific impulse of TH 2 /LOX versus mixture ratio 
for different area ratios and chamber pressures. 


The characteristic velocity is shown in Figure 4. For the approximation 
of c* the expression c*(l + l/r) can be used, which shows in the near of the 
optimum a fairly linear characteristic [3]. Selecting a power law formulation 
for the chamber pressure influence [4] , the subsequent relation for the charac- 
teristic velocity within the limits 4^ r^ 7 results. 


,* = 


3660 - 160 r 
1 + l/r 


70/ 


- 0. 022 


(m/sec) 


(31) 
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m/sec (ft/sec) 


2500 



2800 N/cm 2 
(4000 psia) 

700 N/cm 2 
(1000 psia) 


Figure 4. Characteristic velocity of LH 2 /LOX versus mixture ratio for 

p = 700 N/cm 2 (1000 psia) and 2800 N/cm 2 (4000 psia) . 
c 


The equation is valid for chamber pressures between 500 N/cm 2 and 3000 
N/cm 2 ; however, both limits may be extended without decreasing the accuracy 
significantly. 

The thrust coefficient is shown in Figure 5. To simplify the approxi- 
mation with the common equation for C a y can be defined, which is 

F o C F 

easier to use, 
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Figure 5. Thrust coefficient of LH 2 /IjOX versus area ratio 
for different mixture ratios and chamber pressures. 


This y is plotted in Figure 6. For space flight applications, the area ratio 
C F 

(A /A ) is preferred over the pressure ratio (p /p ), since the latter value is 
e t c e 

more characteristic for the engine size. The relation for y is 


F 


y 


c 


F 


exp 


0.00534 



- 0. 0311 (r - 3) 



-0. 0555 


(33) 
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A e /A t 


Figure 6. Specific heat ratio y for the calculation 

C F 

of the thrust coefficient as a function of area ratio for 
different mixture ratios and chamber pressures. 

The applicable pressure and mixture range are the same as for the character- 
istic velocity; in addition, the area ratio ranges from 50 to 500, where again 
the limits may be extended. 

Finally the pressure ratio as a function of the area ratio is of interest. 
This correlation is shown in Figure 7. According to Reference 5, a power 
relation of the following type represents the characteristics of both quantities 
well: 



exp 


(1.38 - 5.68 10" 4 r) in 



+ 1.58 


- 0. l(r - 3) 




(34) 
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